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Abstract 

Lock step walker model is a one-dimensional integer lattice walker model in discrete time. Suppose that 
initially there are infinitely many walkers on the non-negative even integer sites. At each tick of time, each 
walker moves either to its left or to its right with equal probability. The only constraint is that no two 
walkers can occupy the same site at the same time, ft is proved that in the large time limit, a certain 
conditional probability of the displacement of the leftmost walker is identical to the limiting distribution 
of the properly scaled largest eigenvalue of a random GOE matrix (GOE Tracy- Widom distribution). The 
proof is based on the bijection between path configurations and semistandard Young tableaux established 
recently by Guttmann, Owczarek and Viennot. Statistics of semistandard Young tableaux is analyzed 
using the Hankel determinant expression for the probability from the work of Rains and the author. The 
asymptotics of the Hankel determinant is obtained by applying the Deift-Zhou steepest-descent method to 
the Riemann-Hilbert problem for the related orthogonal polynomials. 



1 Introduction 

In Jl^ , two types of random vicious walkers models, random turn walker model and lock step walker model, 
are considered. In these models, walkers are on one-dimensional integer lattice, and time is discrete. For their 
applications and earlier results, see, for example, |l5[ |l^, |8[ 0, ^, and references therein. In this 
paper, we present resuhs on lock step model showing a relation to random matrix theory. For random turn 
walker model, see || and discussions following Theorem LI below. 



At time t = 0, infinitely many walkers are located at the sites {0,2,4,6,- • •}. We label the walkers by 
Pi, P2, P3, ■ ■ ■ from the left to the right. In the lock step model, at each time t = n, all the particles move either 
to their right or to their left with equal probability. The only constraint is that no two particles can occupy 
the same site at the same time. This is why the walkers are called "vicious" . One typical path configuration is 
shown in Figure |^. 

This model can also be thought of as a certain totally asymmetric exclusion process in discrete time. 
Initially there are infinitely many particles at {1,2,3,- • •}. A particle is called left-movable if its left site is 
vacant. Particles Pj+i, Pj+2, • • • ,Pk are called successors of a particle Pj at a certain time if they are next to 
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Figure 1: vicious lock step walkers 



each other in the order of the indices. At each (discrete) time step, a left-movable particle either moves to its 
left site together with arbitrarily taken number of its successors, or stays at the same site with equal probability. 
It is easy to see that this process is equivalent to the above lock step model ; right move of lock step corresponds 
staying at the same site in the exclusion process. Figure represent an exclusion process equivalent to the lock 
step model in Figure 
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Figure 2: exclusion process 



Suppose that during N time steps, total k left moves are made by all the particles. In the example of Figure 
|l|, = 6 and k = 14. We denote by P(A'', k) the set of all path configurations during N time steps with total k 
left moves. Then each configuration has equal probability, 1 over the cardinal of P(A^, k). Hence our probability 
space is P{N, k) with uniform probability given by 1/| P(A^, fc)|. We denote by Lj{N, fc)(7r) the number of left 
moves made by the particle Pj in a path tt e P(A^, k). We are interested in the limiting statistics of the random 
variables Lj{N, k) a.s N,k ~^ 00. 
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To state the main result, we need a definition. Let u{x) be the solution of the differential equation 

Uxx — "^u^ + xu, u{x) ^ —Ai{x) as x +00, (1-1) 
where Ai is the Airy function. The above equation is called Painleve II equation. It is known that there is unique 



global solution to (1.1) (see, e.g. B and references in it). Define the function, called the GOE Tracy-Widom 



distribution function, by 

Fi(x) =exp{ - -y (s~x){u{s)fds + -J u{s)ds}. (1.2) 
This is indeed a distribution function, and the decay rate is given by 

Fi{x) = l + 0(e^^^'''), asx^+oo, (1.3) 
Fi{x) = 0(e-"l^l'), asx->-c5o, (1.4) 

for some c > (see, e.g. (2.11)-(2.14) of [Q). In |29), Tracy and Widom proved that Fi is the limiting 
distribution of the properly centered and scaled largest eigenvalue of a random matrix taken from a Gaussian 
orthogonal ensemble. The subscription 1 in Fi is a general convention : there are also GUE and GSE Tracy- 
Widom distribution functions F2 and F4 ||2^ , One can find general discussion for random matrices in 

Now the main theorem is 



Theorem 1.1. For fixed <t <1, let 



,1/3 



Vit) = ^. p(i)^lfc^. (1.5) 
1 + t l + t ^ ' 

Let Fi{x) be the GOE Tracy-Widom distribution function. Under the condition that in N time steps total k 
left moves are made, the (conditional) probability distribution of the number Li{N,k) of left moves made by the 
leftmost particle satisfies 

pi^^l^^l2:iL^!^<^]^F.(x). when k ^ 



rl^A^'^^im^-')^^'^'^^ ft = ^iV^ + 0(iVV^). (1.6) 



Also all the moments of the scaled random variable converge to the corresponding moments of the limiting 
random variable. 

In other words, in the large N limit, the (conditional) fluctuation of the displacement of the leftmost particle 
in lock step model is identical to the fluctuation of the largest eigenvalue of a random GOE matrix. Naturally 
we expect that the k^^ particle corresponds to the /c*^ eigenvalue of random GOE matrix. 

It is interesting to compare the above result with the results for random turn walker model. Initially there 
are infinitely many walkers Qi, Q2, Q3, • • • at the position {1, 2, 3, • • • }. We again call a walker left-movable if 
its left site is vacant. At each time, we select one walker at random among left-movable walkers, and move it 
to its left site. Hence there is one and only one movement at each time and all the movements are to the left. 
An example of random turn walker path configuration is in Figure ^ Let Xj{N) be the displacement of the 
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Figure 3: random turn walker model 



j^^ walker after N time step. It is shown in ||^ that for j = 1, 2, we have 



lim I 



Xj{N) - 2VN 



<x]^Fi'\x), J = 1,2, 



(1.7) 



where Fi'''' is the limiting distribution of the (scaled) j*'' largest eigenvalue of a random GOE matrix. Especially 
we have Fi{x) = Fi^\x). On the contrary, if we assume that the walkers move to their left in the first N time 
steps, and then move to their right in the next N time steps so that at the end walkers come back to their 
original positions, then we obtain the GUE Tracy- Widom distribution in the limit |Q. Indeed in this case, a lot 
more are known. The general j^^ row statistics and also the correlation functions converge to the corresponding 
quantities of random GUE matrix in the limit . 

The first step to prove the above theorem is to map the path statistics into tableaux statistics following . 
By definition, a semistandard Young tableau (SSYT) is an array of positive integers top and left adjusted as in 
Figure ^ so that the numbers in each row increase weakly and the numbers in each column increase strictly. A 
reference for tableaux is p6t , and we freely use the notations in it. In JTsf , Guttmann, Owczarek and Viennot 
established a simple bijection between path configurations of lock step model and the set of SSYT : for a path 
configuration, we write down the time steps at which the j"^ particle made movement to its left on the j*"^ 
column. Hence the top row is the array of time steps the particles made first movement to their left, the 
second row is the array of time steps the particles made their second movement to their left, and so on. If we 
draw boxes around each number, the result is a SSYT. See figure ^ for the tableau corresponding to the path 
configuration of Figure 0. This map is a bijection between P(A^, fc) and the set of SSYT of size k with fillings 
taken from {1, 2, • • • , N}. Moreover, under this bijection, Lj{N, k) is equal to the number of boxes in the j^^ 
column of the corresponding SSYT. Therefore the statistics of Lj is identical to the j"^ column statistics of 
random tableaux. 
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Figure 4: semistandard Young tableau 



After the work of Guttmann, Owczarek and Vicnnot, Forrester in observed that similar bijection can be 
established between path configuration of random turn model and the set of standard Young tableaux (SYT). 
The above result (1.7) is obtained based on this bijection and the recent work M on statistics of the first row 



of random SYT. Also as mentioned above, if we assume that the walkers move to their left in the first N time 
steps, and then move to their right in the next N time steps so that at the end walkers come back to their 
original positions, the limiting fluctuation is not GOE type, but GUE type. This difference comes from the fact 
that in this case, the path configuration is in bijection with the pairs of SYT. The statistics of pairs of SYT is 
more well studied than that of single SYT |2^, ||, 19 , and we have stronger results mentioned earlier. 



This paper is organized as follows. In Section]^, using the result of we express the generating function 
for the probability of the first column of random tableaux in terms of a Hankel determinant. It is a general fact 
that Hankel determinant is related to orthogonal polynomials on the unit circle. The asymptotics of orthogonal 
polynomials is obtained via Riemann-Hilbert problem and summarized separately in Section ^. We can obtain 
the limiting statistics of the first column from the knowledge of Hankel determinant asymptotics. This work 
occupies the second half of Section ^j. The proof of Theorem 1.1 is given at the end of Section |^. 

Acknowledgments. The author would like to thank Percy Deift for his interest and encouragement. 



2 Proof 

Let d\{N) be the number of semistandard Young tableaux of shape A with fillings taken from {1.2. •• • ,N}, 
and let ^(A) be the number of rows of A (parts of A, or the length of the first column). From the bijection of 
path configurations and tableaux, the number of path tt e P{N, k) satisfying Li{N, fc)(7r) < / is equal to 



(2.1) 



Ahfc 

e(\)<i 



In our analysis (see also Q), it turns out that in addition to the number of rows, the number of odd columns 
plays an important role in describing a tableau. For a partition A, we define A' to be the transpose of A, /(A) 
to be the number of odd row in A, and |A| to be the size of A. Let b{N^ j,m,l) be the number of semistandard 
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Young tableaux of size 2j + m with m odd columns with at most I columns with fillings taken from {1, 2, • • • , N} 



biN,j,mJ)^ ^ dxiN). (2.2) 

\h2j+m 
/(A')=m 

e(x)<i 

We use the notation b{N,j, m, oo) for the sum above without restriction on ^(A). Now we define a generating 
function 

i(X)<l 

where the sum in the first expression is taken over all the partitions A satisfying ^(A) < /. 
The starting point of our analysis is the following result of H . 



(2.3) 



Lemma 2.1. Let (j){N,l,t, (3) he defined as in (2.3). We have 

cj){N, 21 + 1, t, /3) = (1 - f2)A^(JV-i)/2 det(i/0, (2.4) 
where Hi = {hjk)o<j^k<i is the I x I Hankel determinant with 

oj+fc+l /•! 

hjk^hjk{N,t)^ / x^+^il + t' -2tx)-^{l-x'f/^dx. (2.5) 

J-i 



Remark. Note that the right hand side of (2.4) does not depend on (3. 

Proof. This proof is in ||^ in a slightly different form. Let sx{x)^ x = (xi,a;2, • • • )i be the Schur function, and 
define H{u] y) with y = (yi, y2, • ■ • ) by 

H{u-y)=\{{l-uy,)-\ (2.6) 

j 

In (5.65) of it is proved that 

P^^^">sx{x)^HiP;x)Eu^spi2i)detiHiU;x)), (2.7) 

e{X)<2i+i 

which is an identity as a formal power series in x. But the combinatorial definition of the Schur function is 
(see, e.g. Chapter 7.10 of p§) 



sx 



{x)^Y.-T'"'-T'^'---^ (2-8) 



where the sum is over all semistandard Young tableaux T of shape A, and aj{T) is the number of parts of T 
equal to j (type of T). Since aj{T) = |A|, if we take the special case x — {t,t, - ■ ■ , i, 0, 0, • • • ) where the first 
N elements are t and the rest are 0, then sx{x) becomes 

t^^^dxiN). (2.9) 
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Hence for this special choice of x (EJ) is now 



^ /3/(^')tl^ldA(7V) = (l-/3tr^Ec,eSp(2odet((l-i[/r^). (2.10) 

£(A)<2i+l 



Now using Weyl's integration formula for Sp{2l) (see, e.g. [g5|), the expectation in (2.10) becomes 

Ec/6Sp(20 det((l - ft/)-") = / TT {cose^~cosekf\\siu'6^{l + e-2tcos6,r^de,. 

(2.11) 

Standard Vandermonde determinant manipulations yield 

2'' / /■^'^ \ 
- — -ydet / cosJ+'= 61 sin2 61(1 + i2-2t cos 6i)-^d6l) , (2.12) 

(271') \7o 



0<j,k<l 



which again after change of variables x — cos 6, is equal to 

2' 



\dei( [ Uj{x)Uk(x){\+t^ -^txY^il-x^f^dx] (2.13) 

\J-l J Q<j,k<l 

where Uj{x) — ^^^^^^^^j^, x ~ cos 9, is the Chebyshev polynomial of the second kind. Note that Uj{x) = 
2^ x^ + • ■ ■ . Hence elementary row and column operations yield Lemma □ 

Using this expression, we first obtain the asymptotic result for the generating function. The limit is insen- 
sitive to (3 since so is cj). 

Proposition 2.2. Let < t < 1 and f3 > Q be fixed satisfying < (3t < 1. For each I and N, define x ^M. by 

l~r]{t)N 



p{t)m/3 



(2.14) 



where ri{t) and p{t) are defined in (L5) Then there exits a positive constant Mq such that for M > Mq, there 
are constants C, c > 0, independent of M , and C{M) which may depend on M so that 

\cj,{NJ,t,(3)-Fiix)\<^^ + Ce-'^'''\ -M < x < M. (2.15) 



Also we have 



< 1 - 0(iV,/,i,/3) < Ce"^"^''", x>Mo, (2.16) 

< 0(iV, I, t, P) < Ce-'^l'^l', X < -Mq. (2.17) 



Proof. It is enough to consider the limit for <j){N,2l + l,t,P) since from the definition (2.3), is monotone 
increasing in I. First we related the determinant in (2^) with certain quantities of orthogonal polynomials on 
the circle. 

Let Pj{x) — x^ + ■ ■ ■ be the j^^ monic orthogonal polynomial with respect to the weight w{x)dx = (1 + — 
2te)~^(l — x^y^^dx on the interval (—1, 1), and let Cj be the norm of pj : 

r-l 

P]{x)pk{x)w{x)dx = CjSjk- (2-18) 
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It is a well known result of orthogonal polynomial theory (see, e.g. |2j) that Cj = dct{Hj+i)/ dct{Hj), where 
Hi = {hjk)a<j,k<i with 

hjk = / x^^'^il + f- 2tx)-^{\ - x^f/^dx. (2.19) 



which is equal to J^+i hjk (recall ( |2.5| )). Hence det(iJj) = ^det(iJj). Since the Szego strong limit theorem 
for Hankel determinants (see, e.g. [^) implies that limj^oo det(_ff;) = (1 — ^2^-^(^-1)72 fixed N, we have 

2, + 1, W) = 1™, fl = n . (2.20) 

j = l J T / ^-^J 

Now we use the relation between orthogonal polynomials on the unit circle and those on the interval (—1, 1). 
Let ■nj[z) = + ■ ■ ■ be the j*"^ monic orthogonal polynomials on the unit circle {j^l — 1} with respect to the 
weight 

^{z}^ = (1 - tz)-^{l - tz-^''^, (2.21) 



and let Nj be the norm of 7Tj{z) 



dz 



^jizWizMz)—- = NjS.k- (2.22) 
\z\=i 27riz 



There is a simple relation between orthogonal polynomials pj on the unit circle and orthogonal polynomials tTj 
on the interval (see the forth equation of (11.5.2) in [^) : 

Cj ' Pj(x) = W-(l + 7r2j+2(0)) iVa^Yi , , x ^ -{z + z ). (2.23) 



TT Z ^ Z 

Especially comparing the coefficient of the leading term , we have the relation 



= + ^2j+2(0))7V2j+i. (2.24) 



But we also have (1 - 7r„+i(0)2)iV„ = A^„+i (see (11.3.6) in ||27|1). Hence ( |2.24| ) is equal to 

71 

227 



Cj^l^{l-7T2j+2iO)) 'N2J+2. (2.25) 



Therefore (2.20) becomes 



</)(7V, 2/ + 1, t, (3) = - vr2,+2(0)) A^2-^2- (2-26) 



This argument appeared in Corollary 2.7 of Now using Proposition 3.2, computations similar to Lemma 
7.1 of B (also Corollary 7.2 and Corollary 7.6 of [i) yield the result. □ 



Interpreting the notation 6(iV, j, m, cx)) as the sum without constraints on 1{X) in (2.2), the number of path 
configuration in P(iV, fc) is equal to 

|P(7V,fc)|= HNJ,m,^), (2.27) 

2j-\-m—k 
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and the probability of interest that Li{N, k) < I for tt e P{N, k) is equal to 



|P(iV,fc)| 



2j+m— fc 



2j-t-m— A; 



where 



p{NJ,mJ) 



b{N,j,mJ) 



(2.28) 



(2.29) 



6(7V, j, m,oo) 

For fixed N, sls I ^ oo, the Szego strong limit theorem for Hankel determinants (see, e.g. |2Q]) implies that 



(vM becomes 1. Thus we have the identity 



J2 t^'iPtrb{N,j,m,^) = il-t'y 

j,m>0 



7VC7V-l)/2 



(l-/3<) 



-N 



(2.30) 



By taking Taylor expansion of the right hand side in t and we obtain 

'N{N-1) 



b{N, j, m, oo) 



+ j-l\(N + m-l 



J 



(2.31) 



There is a more direct way to see this. See the remark after Lemma 2.6 below. Now from (2.27), the total 
number of paths in P(A'', k) is equal to 



2j+m— 



— ^^"2 — - +j — l\jN+m~l 
3 



(2.32) 



Now it is straightforward to obtain the following result on the number of all paths. 
Lemma 2.3. Let < t < 1 and let 



As N oo, we have 



^tN{l - t)^(l - f2)JV(JV-l)/2-l 

where the term o(l) vanishes as N ^ oo, and fi is defined by 

k /2 

N 



N l-t2 



(2.33) 



(2.34) 



(2.35) 



which is of order o{N^/^) from (2.33). Moreover, the main contribution to the sum comes from \m — j^N] < 

^l/2+./2 

some < e < i ; precisely, there is a constant c > such that for any < e < we have 



P(iV,fc)| = 



k—m is even 



(l + 0(e-^^ )). 



(2.36) 
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Proof. From (2.32), we have 

[#] 

|P(iV,fc)| = a{m), a(m) 
The ratio of a{m) is 



■m=0 



2 ^ — l\fN + m-l 

/ V m 



a{m + 2) {N + m + 1){N + m){k - m) 

a{m) ~ [m + 2)(m + 1){N{N - 1) + fc - m - 2) ' 



(2.37) 



(2.38) 



One can directly check that under the condition (2.33), the above ratio is decreasing in m, and becomes closest 
to 1 at 



t 



(2.39) 



Hence a(rn) is unimodal: it is increasing for m < rric and is decreasing for m > m^. Now consider the 
neighborhood J\f of nic of size iV^/^+'^/^ for some fixed < e < i. For m in Af, set 



1-t 



N + x, \x\<N^^^+'/\ 



For any M, a; > 0, Stirling's formula yields 



{M + x)l = V2^M''+-e~''+& I 1 + O(^) + O(^) 



Using (|2.33D , ( ^lo| ) and ( ^.4l[ ), we have for m in TV, 

^7V + m-l\ ij^^N + x-iy. exp{-(^ + :E)logi- iig]^} 



(2.40) 



(2.41) 



m J {N ~ ly.ij^N + xy. 



l + 0(A^-2 + 2^) (2.42) 



and 



N(N-l) . k- 



2 ' 2 

fc— m 



(f^ + (M-TjF)f-f-l)! 
(^-l)K^l^ + (M-T^)f-f)! 



Thus we have 



a(m) 



exp{-(^jV^ + (^ - ^)iV - X) \ogt- - ir} , 



exp{-a^ logt - iV^logt - - 1)^ - ^} 

(1 + o(l)). 



(2.43) 



Let 



\/27rt3/2iV3/2(l _ i)JV(l _ t2)Ar(Ar-l)/2- 

|P(7V, fc)H^a(m) + ^a(m), 



(2.44) 



(2.45) 



where * denotes the set Af of m satisfying (2.40) and denotes the rest of the range of m. From (2.44), the 



first sum over * is equal to the right hand side of (2.34). Also from the unimodality, a(m) in is less than or 
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equal to the largest of a(TO+) and a(m_) where m± = [j^iV ± iV^/^+'^]. The number of summand in is of 



order N^. Hence using ( ^.44 ) again for m±, if we take 



4t ' 

-cN 



for large N, we obtain 



which establishes the proof. 



Now we rewrite (2.3) as 



^{N, l,t, t2)^(A'-i)/2 (1 _ f^t)N ^ ^2, ^fitrb^N, J, m, ^)piN, J, m, I). 

The asymptotics of (piN, I, t, (3) and p{N,j, m, I) are related as follows. 

Lemma 2.4. For any d > 0, there are constants Co, Cq > such that for all I > 0, 

p{N, /x+(iV), ,y+{N),l) - < (/.(TV, I, t, P) < p{N, /x_ {N),,y_{N),l) + ^ 



(2.46) 
□ 

(2.47) 



whe 



y±{N) 



t^ 



2(1 -t2) 
I3t 



1 -/3t 



TVicoV^VlogiV 



(2.48) 

(2.49) 
(2.50) 



The proof follows by using the following Lemma twice for j and m indices together with the Lemma 2.6 



(Recall the (|2.3lD). 

Lemma 2.5. For a sequence {Qj}j>o, we define the following generating function 



GiN) = il-a)^J2' 



3=0 



N + j-l 

j 



9j> 



< a < 1, iV = l,2,--- 



(2.51) 



For each d > 0, there are constants Ci,ci > such that for any sequence {<Zj}j>o satisfying (i) q.j > qj+i anc 
(a) 0<qj < 1, 



whe 



N* = 

N** = 



1 - a 
a 



N ~Ciy/N\ogN, 



N + ci^/N \ogN. 



Proof. This proof is parallel to that of the de-Poissonization lemma in pl| . We have 



G(iV) /, = (l-a)^a^Y^+:' ^ 



(2.52) 

(2.53) 
(2.54) 

(2.55) 
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Stirling's formula yields for n,m> 1, 

<Cexp<^m 1 + — log 1 + - log- 

\ m J [_ m m m n 

with some constant C. Thus we have 

fj < Cexp{Nh{j/N)}, h{x) = {I + x) log(l + x) - xlogx + xloga + log(l - a). 
One can directly check the following estimates of h : 



Q<x < 



2a 1 — a 

h{x) < -[log2 - i-i-^ log(l + a)]a;, x> 
2a 



2a 
1^' 



2a 
1 - a' 



We take a constant ci > satisfying — 1 > d. From ( 2.58 ) and the condition (ii), we have 



with a new constant C. Similarly, 



N"<j<-^N 



Also using ( p.59| ), we have 



for some constants c', C". Thus we have 



m) - E f^i^ 

N'<j<N'* 



< 



with a possibly different constant C. Now from the monotonicity condition (i), we have 

E fjlJ^i E fjjiN' <1N', 



N'<j<N* 



N'<j<N' 



and 



E ^ ( E ]<1N" > (1 - T7d)9iV" > qN" - -TT^, 



N'<i<N" ^N'<j<N-" 

using the equality ( 2.63| ) for the second equality. Thus we obtained the desired result. 



(2.56) 

(2.57) 

(2.58) 
(2.59) 

(2.60) 

(2.61) 

(2.62) 

(2.63) 
(2.64) 

(2.65) 
□ 



To use the above Lemma to (p, we need monotonicity in I. It is more convenient now to view semistandard 
Young tableaux (SSYT) as generalized permutations. A two-rowed array 



k3i 



Ik 



(2.66) 
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is called a generalized permutation if cither < v+i or ij. = jr < jr+i- Suppose the elements in the 

upper row of tt come from {1, 2, • • • , M} and the elements in the bottom row come from {1, 2, • • • , N}. One 
can represent a generalized permutation as a M x iV matrix {aik) where aik is the number of times when (^) 
occurs in tt. For example, the generalized permutation 




(2.67) 



corresponds to 



I 





2 Q\ 





3 


1 







1 V 



(2.68) 



In the proof of the Lemma below, we regard a generalized permutation as a M x iV square board with stacks 
of balls in each position (i, k). We denote by L{n) the length of the longest strictly decreasing subsequence 



of TT. In the example (2.67), L{'k) — 2. 

Let Mj^m be the set oi N x N matrices tt 



{ttik) which is symmetric aik = olm, and satisfies ^ 



N 



= m 

and X]i<i<fc<ArOifc = j- This is a certain subset of the set of generalized permutations. The celebrated 
Robinson-Schensted-Knuth correspondence |^ establishes a bijection between Mj^m and the set of SSYT of 
size 2j + m with m odd columns with fillings taken from {1, 2, • • • , iV}. Moreover, under this bijection, L{'k) 
for TT e Mj^m (viewed as a generalized permutation) is equal to the number of rows of the corresponding SSYT. 
With this preliminary, we can prove the following. 

Lemma 2.6 (monotonicity). For any j,m > 0, we have 

p{N,j + l,m,l)<p{N,j,m,l), p{N, j,m + 1,1) < p{N, j,m,l). 

Proof. We first consider the second inequality. From ( |2.29| ), we need to show that 

(to + l)b{N,j,m + 1,1) < {N + m)b{N,j,m,l). 



(2.69) 



(2.70) 



By the definition (2.2) and the Robinson-Schensted-Knuth correspondence, b{N,j,m,l) is equal to the number 
of TT e Mj^rn Satisfying L{ti) < I. 

Consider all possible distinct (strict) upper triangular parts of elements in Mj^m- It is equal to putting j 
identical balls into N{N — l)/2 boxes ; K = (^^(^^i)/2+i-i~j distinct ways. Hence we have a disjoint union 
Mj,m — ^iLiSi,m where each Si^m consists of tt e Mj^m with same upper triangular part, and elements in Si^m 
and Sii ,m have different upper triangular parts when i ^ i' . Similarly, Mj^m+i = ^'^iS'i^rn+i where a G S'i^^y^j^i 
has the upper triangular part same as that of tt e Si^m- 

Now for each tt = (a^s) G 5'i,,„, we generate N + m elements in S'i „^_f_i as follows. For 1 < r < TV, assign 



+ 1 identical tt' = (a'^;) such that a[ 



1 and a'f.[ = aki for {k,l) ^ {r,r). (One can think this as 

1 = m -I- A'', there 



adding a new ball in an array of a^r balls ; there are a^r + 1 ways.) Since J2i<r<N ^rr 
result (to -|- N)\Si^rn\ (many identical) elements of Si_ra+i- Note that under this assignment, 

L{ti') > L{-k). 



(2.71) 
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Now fix cr = (bki) G 5*1, m+i- Since there are m + 1 diagonal entries, there are exactly m + 1 (many identical) 
elements of Si^m from which a is generated under the above assignment. (Considering each entry as a ball, each 
m+1 balls on the diagonal can be a newly added one.) Thus we have the identity (■m+N)\Si^m\ = i'n'i+l)\Si_m+i\- 
Furthermore, from the remark regarding ( 2.7l| ), we have (m + l)|i?(«, m + 1, 01 ^ + N)\Ri^,n^i\, where Ri.,n,i 
is the subset of tt S S'i,m satisfying L(tt) < I. Therefore the second inequality in the Lemma is obtained. 

The first inequality follows from a similar argument. □ 

Remark. As mentioned before, using the generalized permutation interpretation of SSYT, we can see ( ^.31 ) 
directly. From non-negative integer matrix representation of generalized permutations, b{N, j,m,oo) is the 



number of N x N matrices (a^s) with non-negative integer entries such that '^rr 
j. It is equivalent to placing m identical balls of color 1 into N boxes and j identical balls of color 2 into 
N{N - l)/2 boxes. Therefore we obtain ( |2.3lj ). 

Now we give the proof of the theorem. 



and J^Kr 



<s<N ^rs 



proof of Theorem 1.1. In ( ^.28 ), we have 
P(Li(7V,fc) < = 



1 



P{N,k)\ 



^ p{N,j,m,l)b{N,j, 111,00). 



(2.72) 



2j+7n—k 



where P{N,k), p{N, j,m,l) and b{N,j,m, 00) are given in ( 2.32 ), ( ^.3l|) and ( ^.29|) , and b{N, j,m,l) is given in 
( ^.2| ). We split the above sum into two pieces. One part is the sum over (1) |m — jzrt^l ^ N^l'^+'^l'^ and (2) 
the rest, where < e < i is fixed. Then since < p{N, j,m,l) < 1, we have from ( 2.36| ) 



\P{N,k)\ 



PiN, J, m, l)b{N, J, m,^) = 0(e-^^' ) 



(2.73) 



(2) 



for some c > 0. 



We use Lemma 2.4 to estimate p{N, j,m,l) for {j,m) in (1). Set 
k-m- 2coN^logN 



P 



+ k-m- 2coN^/]ogN 



f3 



m — cq^/N logN 



N + m- Coy/ N log N 



where k satisfies the condition in (1.6), and we take t > 0. For (j, m) in (1), they satisfy 



(2.74) 



(2.75) 



The first inequality of ( |2.48| ) yields p{N, j, to, I) < (j){N, + CqN-'^. Set x by (|44|) where t is replaced by i 

and I is given by I = [r]{t)N + xp{t)N^/^]. Let MIMq satisfies — M < 2x < M where Mq is given in Proposition 
23. From (|2.75D, we have 



X = X + 0(1), 



Proposition 2.2 implies that 



|0(iV,Z,i,/3)-Fi(£)| < 



C{M) 
F/3 



I = [r]{t)N + xp{t)N^^^] 



(2.76) 



(2.77) 
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for large N. Since F{{x) — —^{u{x) + v{x))Fi{x) is bounded for a; G M, we have from (2.76) that ^1(0;) = 
Fi{x) +0(1). Thus we have for large N, 



piN,j,m,l) < +0(1), 



(2.78) 



where o(l) term is independent of {j,m) in (1) and vanishes as — s- 00. Thus using Lemma 2^, we have for 
large N, 



(iV, fc) < < Jp^^\ E(^i + <^MN, J, m, 00) + 0(e-^^" ) 



(1) 



= Fi(x)+o(l), l^[r,it)N + xp{t)N^^^] 
Similarly, we obtain the lower bound using the second inequality of ( 2.48| ). Thus we proved ( |l.6| ). 



The convergence of moments is also similar using (2.16) and (2.17) (cf. Section 8 of 



(2.79) 



□ 



3 Asymptotics of orthogonal polynomials 



This section is devoted to asymptotics of orthogonal polynomials used in the proof of Proposition 2.2. The key 
ingredient is the equilibrium measure (see jlO, 11, ^). 

On the unit circle, the equilibrium measure d^v{z) = V'(^')|f for V{z) and its support are uniquely deter- 
mined by the following Euler-Lagrange variational conditions : 



there exits a real constant / such that, 

2 / log \z - s\d^.v{s) - V{z) + / = for z e J, 

2 / log \z - s|d^v(s) - V{z) + / < for z e S - J. 

Lemma 3.1. Let 7 > and < i < 1, and let 

V{z) = 7log(l-iz)(l-tz-i). 
Then their equilibrium measure '>p{0)d0/2T: is given as follows. 
• When < 7 < we have J = E, 1^0, and 

7(1-0 



1-7+ 



= 1 + 7^ 



1 + <2 - 2t cos ( 



1-tz 1-tz- 



When 7 > i±*, J = {e*^ : 16*1 < 9^}, where sin^ % 



(l-t)^(27-l) 



4t(7-l) 



0<9c<n, or 



7(1 -t) 
7-1 ' 



(3.1) 



(3.2) 



(3.3) 
(3.4) 



(3.5) 
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and 



(27-l)(l-t) (27-l)(l-t)^ 
I = 27 log — -T log ■ 



and finally 



2(7-1) 

4(7 - 1) COS f 



4i(7- 1)2 



(1-*)" 



4 sin^ f 



2 .2 

sm sm 



(7-l)(l + z) 



(3.6) 
(3.7) 

(.-t)(.-F^v/(-0(-e-)., -e- (3.8) 

where {z — ^)(z — $^"'")+ denotes the limit of z from inside the unit circle. And in this case, the inequality 
of the second variational condition in (3.1) is strict for z G S \ J. 

Proof. The proof given here is similar to the proof of Lemma 4.3 in j^], whose main ingredient is the following 
results of Lemma 4.2 in [Q. Let d^{s) — u{9)d9 be an absolutely continuous probability measure on the unit 
circle E and u{9) = u{—6). Define 

(3.9) 



g{z) = / log(z - s)d^i{s), 

where for fixed s — e*^" £ S, log(z — s) is defined to be analytic in C \ ((— cxd, —1] U {e*^ : —it < 9 < 9o}), and 
log(z — s) ~ log z for z +00 with z e M. Then for z = e*"^ G S, we have 

(3.10) 
(3.11) 

Also evenness of u yields g{Q) — iri. 

• When < 7 < ■ By residue calculation, it is easy to check ipifi)^ = 1- Define g{z) = J^^ log(z — 



9+iz)+9-iz) = 2 \og\z-s\dfi{s)+i{(j) + TT), 
g+{z)-g_{z) = 27ri / u{9)d9. 



e )ip{9)^ as in (3^). Then by direct residue calculation, we obtain 



9'{z) 



-7* 

l-tz ' 



\z\ < 1, 



Since (7(0) = ni and ^(z) '-^ log(z) as z ^ 00, we have 



kl>l. 



(3.12) 



5(2) 



7log(l - tz) + TTi, \z\ < 1, 
7log(l - tz^i) + logz, |2:| > 1. 

Thus, from ( 3.10| ), the variational condition (3T) is satisfied with J = S and / — 0. 



(3.13) 



When 7 > : Set P{z) ~ ^(z — ^)(z — ^ ^) which is analytic in C \ J and /3(z) ~ z as z ^ +00 with 
z e M. Then we have 



/3(0)--l, (3{t)^~\^-t\, p{t-')^-\^-t\. 



(3.14) 
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First, direct residue calculation using ( |3ll ) shows that jl'g^i'iO)^ = 1, hence ip{0)d0 / {2Tr) is a prob- 
ability measure. Now we define g{z) = log(2: — as before. Using (3.14) again, residue 
calculations yield 



7* A _ (7-l)(2 + l) . . 
1-tzy 2z(z-t)(z-t-i)'^^^'^' 



(3.15) 



Thus we have for \z\ > 1, z ^ (— oo, —1) U [t ^, oo), 



g{z) = Jlog(l-te)(l-iz-i) + ilog^ 



and for \z\ < 1 , z ^ (— 1 , , 



g{z) = ^ log(l - tz){l - tz-') + - logz 



^^^^^:^^^mds^g.il)-,^o,il-t), (3.16) 



1 2i(i:')(:-'-^) ^(-)-^- + - - ^'-''^ 



Now we compute .g+(l) + .g-(l). From ( 3.10| ) with z = 1 and the evenness of "0, we have 

8(7-1) 



g+il) + g^il) - TTi 



log 2 sin - 



2 / (1-*)^ 



t ■ 4sin- 2 



,, sin^ - sin^ -d9. 
2 V 2 2 



Substituting x — (sin ^) ^ sin |, the above becomes 
4(7-1) f' 



log 2 sin —X 



2- 



j2 p2 



da;, p 



1 - t 



2Vi 



sm ■ 



7-1 
V27^' 



(3.18) 



(3.19) 



Using — - Jp ^2"^p2 = 1 which is a consequence of the fact that il){6)d6/ {2ti) is a probability measure, 
we have 



g+(l)+.g_(l)-7ri = 21og(2sin^ 



4(7-1) r 



log(a:;) 



rj.2 _j_ p2 



(3.20) 



Now we use /J^ log x ^f^^-i — Jq^^ \og{sin9)d9 = — f log 2 to rewrite the last term in the above as 

4(7-1) (1+^2) /•! logx 



2(7-1) log 2 + 



Also one can show by residue calculations that 



(x^ + p'^)\/l — x-^ 



zdx. 



logo; 



But one can directly verify that 

dx 



dx — ~ — 



dx 



TT logp 



2 71 (x2+p2)^^2— 1 2py/T+p^ 

p^/x'^ - 1 + ^1 +p2 



(x2 - 1 2py/l+p- 



:l0g 



a;2 — 1 — ^Jl+jPa 



C. 



Thus from ( |3.20| ), (|3.2lD and the definition of p in ( |3.19D , we have 

.,,l, + „_(l)-„^2.,„g^ + ,„g<^'-"»-')^ 



27-1 



4i(7- 1)2 



(3.21) 



(3.22) 



(3.23) 



(3.24) 
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which is equal to 27log(l — t) — I. Now for z E J, from (3.16) and (3.17), we have 

9+ (z) + g-{z) ^ V{z) + \ogz ~ I + TTi. 



(3.25) 



Thus (3.10) yields that the first variational condition (3.1) is satisfied. On the other hand, for z e S \ J, 
g+{z) + g-{z) is equal to the right hand side of ( 3.25| ) plus 

where Ci = {e** : < ^' < argz} oriented from ^ to z if argz > 0, and Ci — {e*^ : argz < < — 
oriented from ^^-^ to z if argz < 0. For 6c < argz < tt, ( ^.26| ) is equal to 

5^ V2(7- l)cos! 



t + t-i- 2 cos 
and for —n < argz < —Oc, it is equal to 

'■^ V2(7 - 1) cos 



- \/ cos 9c — cos Odd < 0, 



argz 



t + t-^ - 2 COS 



^ \/ COS 9c - COS 9dd < 0. 



Thus the second variational condition of (O) is satisfied and the inequality is strict. 

For fixed < < < 1, we define a weight function on the unit circle S by 

ip{z; N) := (1 - tz)-^(l - tz-^)-^ . 



(3.27) 



(3.28) 



□ 



(3.29) 



Let 7r„(z; N) = z" + • • • be the n"^ monic orthogonal polynomial with respect to the measure (p{z; N)dz / {2'Kiz) 
on the unit circle, and let Nn{N) be the norm of 7r„(z; N) : 



^„(z; A^)7r„(z; iV)(l - tz)-'^ (1 - tz-') 
We also define 

<(z;7V) = z"^„(z-i;7V). 
Define a 2 x 2 matrix- valued function F of z G C \ S by 



dz 
2'Kiz 



Y 7-n- N\ ■= n\ , J J2 s-z 



n> 1. 



(3.30) 



(3.31) 



(3.32) 



Then Y{- ; n; N) solves the following Riemann-Hilbert problem (RHP) (see Lemma 4.1 in B) : 



Y{z; n; N) is analytic in z G C \ E, 

/l ■^^p(z;N)' 
Y+iz;n;N) = Y-iz;n-N) ^ > 

y(z;n;7V)(^-7^o.^ =/ + 0(i) as 2 



on z €E S, 



(3.33) 
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Here the notation Y+{z]n : N) (resp., Y^) denotes the hmit of Y{z' \n:,N) a& z' z satisfymg \z'\ < 1 (resp., 
\z'\ > 1). Note that n and N play the role of external parameters in the above RHP. One can easily show 
that the solution of the above RHP is unique, hence (3.32) is the unique solution of the above RHP. This RHP 
formulation of orthogonal polynomials on the unit circle is an adaptation of a result of Fokas, Its and Kitaev 
in where orthogonal polynomials on the real line are considered. 



From (3.32), it is easy to check that 

iV„_i(iV)^i = -Y2ii0; n;N), 



TTniz;N) ^ Yniz;n;N), 

n*Jz;N) = z"rn(z-i;n;iV) = y2i(2;n+l;iV)(y2i(0;n + l;iV))-i. 



(3.34) 
(3.35) 
(3.36) 



The asymptotics of the above quantities can be obtained by applying Deift-Zhou method for Riemann-Hilbert 
problem ( [3.33 ). A reference for Deift-Zhou method is In and similar asymptotics are obtained 
for different weight function e*^^+^ ■* as t,n oo. It is interesting to compare the following results with 
Proposition 5.1 in B. 



Proposition 3.2. For each n and N , define x G 

2t N 



by 



1 + t n 



1 



1 - t 
2(l + t)2 



1/3 



Also let 



2N {2N-n){l-t) n{2N ^n){l-tf 

~ °^ 2{N - n) 4i(iV-n)2 



and for — > 

n 



^ - i+* let 



2t 



sm ■ 



(1-t) n{2N-n) 



2 ' 2{N -n) V t 
There exits AIq > such that as n,N ^ oo, the following asymptotic results hold. 



(3.37) 



(3.38) 



(3.39) 



(i)- IfO< TTt^ <aforO<a<l, then 



|iV„_i(7V)- 



II 



|7r„(0;iV)| <Ce^ 



(3.40) 



(a). Ifa<r2L.^<l- Mn^'^/'^ for some M > Mq andO <a<l, then 



for some constants C , c which may depend on a. 

2t N 
1+t n 

|iV„_i(7V)-i - 1 
for some constant C and c depending on M . 
(Hi). If -M < X < M for some M > 0, then 

iV„„i(7V)-i -1- 



k„(0;7V)| < 



C 



(3.41) 



"2(l + t)2- 


1/3 / \ 




1-t 


ni/3 





^„(0;iV) + (-l)^ 



2{l + ty 



1 -t 



1/3 



u{x) 



,1/3 



< 



C 

7m 



(3.42) 



for some constant C depending on M . 
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(iv). If 1 + <^^<a for some M > Mq and a > 1, then 



-A^„_i(iV)-i-l 



(-1)^ 



7r„(0;n;iV)- 1 



< 



G 



(3.43) 



for some constant C depending on M . 



(v). If a < Z*"" some a> 1, 



sm ■ 



-N^-i{N)- 



1 



cos 



7r„(0;n;7V)-l 



< — 

n 



(3.44) 



/or some constant C depending on a. 
Remark. From the calculations analogous to Section 10 of M, in addition to the above asymptotics results, we 



can obtain more results similar to those in Section 5 of For example, suppose x defined in ( 3.37 ) above 
satisfies ci < x < C2 for some constants ci , C2 G K (hence we are in the case (iii) of above proposition) . For 
a > 0, define w G R by 



a 



2(1 +t)- 
1 -t 



1/3 



2w 



Then we have for w > fixed, 

lim (— 1)"(1 + to)^^7r„(— a; A^) = — 11112 {—iw;x), 

n — >oo 

lim (1 + to)^^7r* (— a; A^) = m22{—iw;x), 



and for w < fixed, 



lim (-l)"(l + ta)-^7r„(-a;Ar) = mii(-m; x)et'"'-2^"', 
lim (l + to)-^<(-a;iV) = -m2i(-«w; a;)et"''-2^'". 



(3.45) 

(3.46) 
(3.47) 

(3.48) 
(3.49) 



where m(z; x) is the solution to the Riemann-Hilbert problem for Painleve II equation with special monodromy 
datap = -q=l,r = (see, e.g. (2.15) of 0]). These results are parallel to (5.21), (5.22), (5.25), (5.26) of §. 



But in this paper, we only need Proposition 3.2 above. 

We are not going to present the detail of the proof because the computation is parallel to that of Lemma 5.1 
and Lemma 6.3 of Q| (see also Proposition 5.1 of Q) where the authors consider the same asymptotic problem 
with different weight function 5^'^°''^. Instead we give some indication why we have Painleve II function in 
the case (iii). Let us consider only when a; > 0. Define a 2 x 2 matrix valued function m{z) by 



m{z) 




~{\~tz)-^ 
{\-tz)^ 
z-"(l-tz-i)-^ 

z"(l-tz-i)^ 



\z\ < 1 



\z\ > 1, 



(3.50) 
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where Y is defined in (3.32). From the RHP for Y, m solves a new RHP : m is analytic in z e C \ S, 
m(z) = / + O(-) as z cx), and satisfies the jump condition m+(z) = m_(z)w(z) on z G S where 



v{z) 



z-"(l-tz)^(l-iz-i)- 



■N 



-z"(l-tz)-^(l-tz-i) 




(3.51) 



The choice of m above is related to the equilibrimn measure in Lemma 3.1, The role of equilibrium measure 
in RHP for orthogonal polynomials is discussed in ||l^, |l^ (see also [Q). The two RHP's for Y and m are 
algebraically related and are equivalent in the sense that a solution to one RHP implies a solution to the other 
RHP. Note that the jump matrix v{z) has the factorization u(z) = (z)u_|_(z) where 



1 -z"(l-tz)-^(l-iz-i)^ 



u_(z) 



,z-"(l-tz)^(l-tz-i) 



v+{z) = 



.0 



1 



(3.52) 



Hence by usual deformation technique of RHP, we can bring the matrix v+ to a contour in |z| < 1, and the 
matrix u_ to a contour in |z| > 1. By the assumption of N and n in case (iii), except in a neighborhood of 
z = — 1, we can find a new contour where the off diagonal entries of v± decay exponentially as rt — > oo. Hence 
the main contribution to the RHP as n — > oo comes only from a neighborhood of z = —1. This is exactly 
related to the fact that the support of the equilibrium measure is Lemma 3J. has special point z — —I at which 
a new gap opens up when 7(= ^) = Now let us focus on the neighborhood of z = —1. The (12) entry of 
V is — e'''^^-' where 



h{z) = nlogz - iVlog(l - fz) + iVlog(l - tz"'^). 
1 + s. Using the definition of x in ( ^.37| ), expansion of ( 3.53| ) becomes 



(3.53) 



Set z 

h{z) — nlog(— 1) — X 



1 - t 
2(l + <)2 



1/3 



nlog(-l) - 2xs + -s^ 
o 



(n 



as n 



1/3 < 



2ni/3 



1 - t 



2(1 + ^2 



1/3 



1 - t 
6(l + <)- 



1/3. 



+ ■ 



(3.54) 



where s ^ [ 2(1+*)^ ] {n^^^ s). Thus we are lead to a RHP with the jump matrix 



i{z) 



f-l)"e 



-2(-a 







(3.55) 



on the line oriented from +ioo to ~ioo. After rotation by — 7r/2, this is precisely the jump matrix for the 
Painleve II equation with parameter p = — (7 = 1, ?■ = (see, e.g. (2.15) of Q : the term (—1)" in off diagonal 
entries can be simply conjugated out). Thus the m, therefore 1", can be expressed in terms of Painleve II 
solution in the limit n — > cx) in the case (iii) with x > 0. 
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